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Abstract 

The third-order elastic moduli and pressure derivatives of the second-order elastic constants of novel 
B2-type AIRE (RE=Y, Pr, Nd, Tb, Dy, Ce) intermetallics are presented from first-principles calculations. 
The elastic moduli are obtained from the coefficients of the polynomials from the nonlinear least-squares 
fitting of the energy-strain functions. The calculated second-order elastic constants of AIRE intermetallics 
are consistent with the previous calculations. To judge that our computational accuracy is reasonable, the 
calculated third-order constants of Al are compared with the available experimental data and other theo- 
retical results and found very good agreement. In comparison with the theory of the linear elasticity, the 
third-order effects are very important with the finite strains are lager than approximately 3.5%. Finally, the 
pressure derivative has been discussed. 

PACS: 71.20.Dg, 62.20.-x, 71.15.Nc 
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1 Introduction 

Aluminium alloys are widely used in the automotive, aircraft and aerospace industries and interest in their applications 

is increasing. Rare-earth element is a class of special elements which are known to be the most effective modifier and 

recrystallizator to the Al alloys. Recently, the Al-TM-RE (TM is a transition metal, and RE denotes a rare-earth element) 

systems have received great attention due to their unique mechanical properties, such as high tensile strength, good 

ductility, high corrosion resistance and thermal stability UJ [2] [3] |4] [5] |6) . Recently, a new class of highly ordered and 

ductile intermetallic compounds with the CsCl-type B2 structure in space group Pm3m (such as YAg and YCu) has been 

discovered with guidance from Gschneider [7| and Morris et al.llSll, It is known that the B2 structure is a very important 
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intermetallics, so the B2-A1RE alloys as well as B2-MgRE alloys have been investigated intensely ll9l ITOl fTTl [121 [131 fT4l 
[15][T6). Some experimental information of structures for Al-RE equiatomic compounds can be found in literature ifTTl . 
and for example, the Al-Y, Al-Ce, Al-Nd, Al-Dy, Al-Pr, etc. are available for metastable B2 phase lfl8l [T9l l20l |2T1 1221 . 
More recently, various studies have been undertaken of the bulk, elastic, thermodynamic, and electronic properties for the 
rare-earth-aluminium B2-A1RE intermetallics from first-principles calculations ll9l fT0irTTl[T2l[T3ll . and these results are 
very important to the novel material design and further scientific and technical investigations. 

For technical applications, applying finite deformation to materials is also very important due to the AIRE inter- 
metallics with superior mechanical properties, while the theory of linear elasticity with assuming infinitesimal strain isn't 
proper any longer and instead the theory of nonlinear elasticity is required [23 . 24, 25, 26, 27). In the linear theory of 
elasticity, the second-order elastic constants (SOECs) are sufficient to describe the elastic stress-strain response and wave 
propagation in solids 1 28 1. In nonlinear elastic theory, high-order elastic constants, such as third-order elastic constants 
(TOECs), play an important role as well as SOECs 11271 - TOECs are not only used in describing mechanical phenomena 
when applying large stress or strain, but also can be used to describe the anharmonic properties such as thermal expansion, 
the interaction of acoustic and thermal phonon, changes in acoustic velocities due to elastic strain, etc.ll29l l30l [3D . In 
addition, the elastic properties related high pressure can also be discussed using the SOECs and TOECs directly. Though 
many experiments have been performed to determined SOCEs and high-order elastic constants G5l . to obtain a complete 
set of TOECs is still not a simple task, especially for crystals with low symmetry or with low-yield stress. Recently, 
a simple method using first-principles calculations have been employed to determined TOECs 021 l33l [34l as well as 
high-order elastic constants, and their results show good agreement with experiments. In this approach, the homogeneous 
deformation strain is applied to the system and usually simple deformation modes are used such as uniaxial tension or 
compression, simple or pure shear, and other combinations of homogenous strains. To provide significant information 
with respect to application and design of these novel AIRE alloys, to study further the nonlinear elastic constants such as 
TOECs of these compounds is required. 

2 Computational methods 

Our method for calculating elastic constants mentions finite-strain continuum elasticity theory. For a solid-body with 
a finite deformation, the configuration of a material point in the system after deformation is represented as x' = x'(x), 
where x is the initial configuration at the equilibrium state. The relation between the strained and unstrained crystal is 
constructed by the deformation tensor 

dx'. 
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where i and y'(=l,2,3)represent Cartesian coordinates. The relation between the elastic constants and the strain energy 
density, O, can be written as 



$0?) = Y\ XI Ci lkl r lij r lkl + X C Uklm„Vij 7 lkl 7 lmn + (2) 
ijkl ijklmn 

where rjjj is Lagrangian strain tensor, which is defined as IT271 

VU = \{Yj JikJ i k ~ 5i \ 



k 

The nf/i-order (n > 2) elastic constants was defined by Brugger [26 1 as 

<9"0> 



Cijklmn... - Cijk... - 



T,=0, 



(3) 



(4) 



where /, J, and K are Voigt subscripts; and the Lagrangian strain tensor T] links this notation by 

V = (»7i»i72,'73»i/4» , 75,%). (5) 

Because of B2-A1RE intermetallics with cubic symmetry, there are three independent SOECs (Cn, C\i, C44) and six 
TOECs (Cm, C112, Ci23,Ci44, C155, C456). To calculate the complete SOECs and TOECs, we introduced six Lagrangian 
strain tensors in terms of a single parameter Inserting these strains into Eq. d2l, the elastic energy per unit mass can be 
written as an expansion in the strain parameter 

= ~A 2 f 2 + + 0{^), (6) 

where the coefficients A2 and A3 are combinations of second- and third-order elastic constants of the crystal, respectively. 
The selections of the different deformation modes leads to different strains used in this work, which are labeled as r) a , 
a = A,B, . . . ,F. The applied strains and the corresponding coefficients A2 and A3 are listed in Table 1 . 

In every case for T] a , £ is varied between -0.08 and 0.08 with step 0.008 to obtain accurate TOECs. The maximal 
amplitude of the deformations is reasonable to obtained the TOECs accurately B2l [331 . The elastic constants will be 
obtained from the least-square polynomial fit to the strain-energy relation from first-principles total-energy calculations. 
In order to obtain the unit cell of strained crystal, the deformation tensor /y is applied to the undeformed lattice vectors 
a,, where i is the lattice index. The deformed crystal is obtained then from aj = Yij Jij a j- T° implement the different de- 
formation modes in our calculation, we need to have the deformation tensor /y, which is determined from the Lagrangian 
strain by inverting Eq. Q, 

J U = Sij + TJij - r 2 mkrtkj + \ 2 I^Mj + " ' (7) 
k kl 

For a given 77, in general, J is not unique but this is not a problem since the Lagrange strain brings rotational invariance. 

We carry out first-principles total-energy calculations based on the density functional theory (DFT) level, using 
the Vienna ab initio simulation package (VASP 4.6) developed at the Institut fur Materialphysik of Universitat Wien 



Il35l [36l [37l . The Perdew-Burke-Emzerhof (PBE) [38, 39 1 exchange-correlation functional for the generalized-gradient- 
approximation(GGA) is used. A plane-wave basis set is employed within the framework of the projector augmented 
wave (PAW) method [40, 41 ]. The equilibrium theoretical crystal structures are determined by minimizing the Hellmann- 
Feynman force on the atoms and the stress on the unit cell. The convergence of energy and force are set to 1.0 x 10" 6 eV 
and l.Ox 10~ 4 eV/A, respectively. For the Billouin zone integrals, reciprocal space is represent by Monkhorst-Pack special 
k-point scheme [42|. Since high accuracy is needed to evaluate the TOECs, the £-point mesh size with 25 x 25 x 25 and 
cutoff energy with 600eV are used to calculate SOECs and TOECs for all calculated materials. Our convergence tests 
show that the chosen parameters are sufficient to reach the desired convergence for the total energy as well as the TOECs. 
In Table 2, we give the equilibrium lattice constants for AIRE interme tallies in our calculation, and it is shown that the 
results agree well with the previous calculation IfTOlfTZl and experiment [18, 19 20l l2Tll22l . 

3 Results and discussion 

For completeness, we also list our calculated SOECs as well as the previous results by Tao and Ouyang, et al. JTOl 
for AIRE intermetallics in Table 2. Our results show good agreement with those obtaining from previous first-principles 
calculations. The bulk modulus B - (Cn + 2C\i)j?> in our calculations also provide extremely good agreement with 
the values fitted from the Rose's equation of state IflOl . To benchmark the reliability results of the presented method, 
we have compared our theoretical results for Al, which is a well-studied cubic crystal having fee structure with available 
experimental data [25 43 1 and previous theoretical findings within DFT theory [34]. The calculated SOECs and TOECs of 
Al are listed in Table 3. In the previous calculations, the all third-order elastic constants of Al were obtained [34]. Overall, 
our ab initio results agree well with the previous calculations and the experimental data for both SOECs and TOECs. We 
present the unknown values of TOECs for all calculated AIRE intermetallics in Table 4. The six strain-energies for A1Y, 
including the results of the first principles calculations and the fitted polynomials, are shown Figure [T| It shows that 
the strain-energies with negative strains are always larger than those with positive strains, so the values of TOECs are 
typically negative. It is worth noting that for the studied intermetallics with Lagrangian strains up to 8.0%, including the 
terms up to third-order in energy expansion sufficed to obtain good agreement with our ab initio results. Next we focus on 
examining for which range of strains the third-order effects dominant the properties of solids. In Figure [2] we also show 
the curves of the linear elasticity comparison with the nonlinear elasticity and the DFT results particular deformation t/a 
in A1Y. It is clearly to see that the third-order effects must be considered and it is not sufficient for linear elasticity when 
applied to finite deformations larger than approximately 3.5%, while the nonlinear elasticity must be considered. 
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4 The effective elastic constants under pressure 

In the case of materials under larger hydrostatic pressure, it is useful to describe the nonlinear elastic properties using 
the concept of effective elastic constants Cy(P). Usually, Cij(P) as a function of pressure P can be expanded as Taylor 
series and it is appropriate to consider only linear term in the external hydrostatic pressure 

Cu(P) * Cjj + '^jjP-P = C u + C tj P, (8) 

with pressure derivatives C' u being material parameters. Once the complete SOECs and TOECs are available, the the 
effective elastic constants under high pressure can be discussed directly. Naturally, the values of C' u are determined from 
components of both SOECs and TOECs, and can be expressed as 1231 

Cm + 2Cn2 + 2C U + 2C\2 



C'u 



2Cn2 + C123 - Cu - C12 



c = - 

C\\ + 2Ci2 

r , _ _2Ci55 + C144 + C\\ + 2Ci2 + C44 

C44= C I1+ 2C 12 ' W 

while the pressure derivative of the bulk modulus is expressed as 

B' = C " +2C ' 2 . (10) 
3 

We give the values of pressure derivatives C' u and B' on the basis of our prediction for second- and third-order elastic 
constants in Table 5. We also list the results of B' obtained from the Rose's equation of state [10], which show good 
agreement with the predicted results from TOECs in our calculations. In Ref.[44|, the following method has been used to 
determine of the pressure dependence of the second-order elastic constants. First, the hydrostatic strain has been applied 
to the crystal, and then the crystal has been additionally deformed to determine the pressure dependent elastic constants. 
The ab initio results for the total elastic energy combined with the strain-energy relation have enabled us to determine 
Cij(P) and C' u . The methods of calculations on the pressure derivatives had also been employed in Ref. ll32l in which the 
related properties have been studied for selected semiconductors. Therefore, we believe that the method used in our paper 
is accurate. Furthermore, the experiments when materials are under high pressure may obtain pressure derivatives, and 
can demonstrate that our results of the nonlinear elasticity for AIRE intermetallics are reliable. 

5 Conclusions 

In this work, we present the unknown TOECs Cjjk and the pressure derivatives C' u of the SOCEs for the novel AIRE 
(RE=Y, Pr, Nd, Tb, Dy, Ce) intermetallics with B2-structure using the first-principles total-energy calculations combined 
with the method of applying a series of homogenous finite strains to the crystals. The PAW method within the GGA 
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is employed. From the nonlinear least-squares fitting, the predictions for the SOECs and TOECs are obtained from the 
coefficients of the fitted polynomials of the energy-strain functions. The calculated second-order elastic constants of AIRE 
intermetallics are consistent with the previous calculations. To benchmark the reliability results of the presented method, 
we have computed the TOECs for Al and compared with the available experimental data and other theoretical results and 
found very good agreement. In comparison with the theory of the linear elasticity, our work shows the nonlinear elastic 
effects must be considered when the finite deformations applied are lager than approximately 3.5%. It is also worth noting 
that for the studied intermetallics and examined range of deformations, including the terms up to third-order in energy 
functions sufficed to obtained good agreement with our calculated results. We have discussed the pressure derivatives 
of the SOECs within the framework of the second- and third-order elastic constants, and the pressure derivatives of the 
bulk modulus agree well with those obtained from Rose's equation of state. We believe that our DFT results of the 
predicted TOECs can be a very useful tool in applying these novel rare-earth intermetallics to practical engineering, in 
which nonlinear effects often really matter. 
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Table 1 The coefficients A 2 and A 3 in Eq .{6]l of corresponding Lagrangian strains as combinations of SOECs and TOECs for crystals 
with cubic symmetry. 



Strain type A 2 A 3 





= (£0,0,0,0,0) 


jCn 




Vb 


= (££0,0,0,0) 


(C„ +C 12 ) 


3C111 + C112 


Vc 


= (£££ 0,0,0) 


|Cn + 3Ci2 


jCni + 3Cn2 + Ci 2 3 


m 


= (£0,0, £0,0) 


jCn + 5C44 


gCnl + ^Cl44 


Vd 


= (£0,0,0,£0) 


\C\\ + JC44 


gCni + 2-C155 


■Hd 


= (0,0,0,£££ 


2C44 


C4S6 



Table 2 Our calculated lattice constants for AIRE (RE=Y, Pr, Nd, Tb, Dy, Ce) intermetallics in comparison with the values of the 
previous calculations and the experiments. The second-order elastic constants (SOECs) and bulk modulus are also listed. 





AIY 


AlPr 


AINd 


AlTb 


AlDy 


AlCe 


a (A) 


3.606" 


3.760° 


3.729" 


3.614" 


3.597" 


3.675" 




3.605*, 3.522 c , 3.754 rf 


3.759*, 3.533'', 3.82 e , 


3.728*, 3.75^ 


3.614* 


3.596*, 3.71* 


3.795*, 3.575 c , 3.85* 


C„ (GPa) 


81.26", 77.93* 


72.61", 66.83* 


74.25", 69.00* 


80.98", 77.13* 


81.40", 78.28* 


62.66", 64.75* 


C12 (GPa) 


54.58", 56.58* 


46.82", 49.26* 


48.71", 51.09* 


55.17", 56.92* 


56.11", 57.73* 


52.64", 47.24* 


C44 (GPa) 


62.79", 61.91* 


47.30", 46.67* 


50.24", 49.83* 


62.30", 61.17* 


64.00", 63.06* 


42.53", 43.31* 


B(GPa) 


63.47", 63.70* 


55.42", 55.11* 


57.22", 57.06* 


63.73", 63.66* 


64.54", 64.58* 


55.98", 53.07* 



"This work 



*ReferencefT0l 
'Reference (T2) 
"■Referenced) 
'Reference |22) 
'Reference] 20] 
s Reference l21l 
''Referenced) 



Table 3. Comparison of the calculated SOECs and TOECs for Al with the previous theoretical values and experimental results. The 
unit of all data is GPa. 





This work 


Previous study 


Experiment 


Cu 


113.21 


110.4" 


108* 


C12 


53.93 


54.5" 


63.92* 


C44 


30.9 


31.3" 


28.3* 


Cm 


-1088.20 


-1253" 


-1080', -1224" -1427"' 


C112 


-415.40 


-426" 


-315", -373',-408"' 


C144 


-14.25 


-12" 


-23', -64',-85"' 


C155 


-479.58 


-493" 


-340", -ses'-sge" 1 


C123 


167.52 


153" 


36", 25',32 £/ 


C456 


-24.18 


-21" 


-30", -27",-42 d 



"Referenced!) 
*Referencel25l 

" At 300K and 298K from References (25) and (43) 
d At 80K from Referenced 
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Table 4. The predicted results of the third-order elastic constants (TOECs) of AIRE (RE=Y, Pr, Nd, Tb, Dy, Ce) intermetallics. The 
unit of all data is GPa. 





A1Y 


AlPr 


AINd 


AlTb 


AlDy 


AlCe 


Cm 


-830.76 


-662.54 


-705.11 


-839.01 


-830.318 


-275.34 


Cm 


-170.00 


-136.39 


-137.27 


-168.03 


-170.36 


-160.98 


C\21 


-200.28 


-202.06 


-231.37 


-199.37 


-210.12 


-264.29 


C144 


-183.22 


-44.48 


-53.57 


-154.28 


-167.32 


97.32 


Cl55 


-355.51 


-290.12 


-300.70 


-348.44 


-352.47 


-269.64 


C456 


-321.66 


-288.12 


-295.93 


-317.28 


-321.05 


-362.77 



Table 5. Predictions for the pressure derivatives of SOECs for AIRE (RE=Y, Pr, Nd, Tb, Dy, Ce) intermetallics. The pressure derivatives 
of the bulk modulus as well as the values obtained from Rose' equation of state 1101 are shown. 





A1Y 


AlPr 


AINd 


AlTb 


AlDy 


AlCe 


C'u 


4.72 


4.19 


4.27 


4.72 


4.62 


2.18 


c 

12 


3.55 


3.57 


3.66 


3.51 


3.56 


4.18 


c 


3.37 


2.47 


2.52 


3.13 


3.17 


1.38 


B' 


3.96 


3.87 


3.89 


3.92 


3.93 


3.57 


B m 


3.99 


3.92 


3.95 


3.92 


3.99 


3.72 
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Lagrangian strain \ 



Figure 1: The strain-energy relations for A1Y. The discrete points denote the values of DFT calculations; solid 
curves represent the results obtained from third-order polynomial fitting. 




-0.08 -0.06 -0.04 -0.02 0.02 0.04 0.06 0.08 



Lagrangian strain % 

Figure 2: Energy as a function of Lagrangian strain parameter £ for particular deformation ^for A1Y. Full 
points denote results of DFT computations; solid and dashed curves indicate the results obtained from nonlinear 
and linear elasticity theory, respectively. 
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